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Abstract The scaling behavior of the thermoelectric trans- 
port properties in disordered systems is studied in the energy 
region near the metal-insulator transition. Using an energy- 
dependent conductivity a obtained experimentally, we ex- 
tend our linear-response-based transport calculations in the 
three-dimensional Anderson model of localization. Taking a 
dynamical scaling exponent z in agreement with predictions 
from scaling theories, we show that the temperature-dependent 
cr, the thermoelectric power S, the thermal conductivity K 
and the Lorenz number Lo obey scaling. 

The scaling description j|] of disordered systems, e.g. the 
Anderson model of localization, has cultivated our un- 
derstanding of transport properties in such systems M-M ■ 
According to the scaling hypothesis, the behavior of the 
d.c. conductivity a near the metal-insulator transition 
(MIT) in the Anderson model can be described by only 
a single scaling variable. As a result of the scaling the- 
ory, the dynamical conductivity in the three-dimensional 
(3D) Anderson model behaves as 

<r(t,T) , / t 



J>l/z Ul/l/Z I ' ' ' 

where T is the temperature and t is the dimension- 
less distance from the critical point. For example, t = 
|1 — Ep/E c \ where Ep and E c are the Fermi energy and 
the mobility edge, respectively. The parameter v is the 
correlation length exponent, which in 3D is equivalent to 
the conductivity exponent, a oc t", and z is the dynam- 
ical exponent, a oc T l ' z . It was further demonstrated 
that not only a(t, T) obeys scaling in the 3D Anderson 
model but also the thermoelectric power S(t,T) |^,0|, 
the thermal conductivity K (t, T) and the Lorenz number 
L (t,T) [Q. However, despite the quality of the scaling 
of cr, we obtained an unphysical value for z J7[|. Scaling 
arguments for noninteracting systems predict z = d in 
d dimensions (|,[5[. But we found | z = 1/v <C 3. In 
addition, values of S(T) (g^ are at least an order of 
magnitude larger than in measurements of doped semi- 
conductors Q and amorphous alloys (ll],[L2j. 

In what follows, we show that we obtain the right 
order of magnitude |u| and good scaling for these ther- 
moelectric transport properties by using a "modified" 
critical behavior of a in the linear-response formulation 
for the Anderson model based on experimental data. 



In the linear-response formulation, the thermoelec- 
tric transport properties can be determined from the ki- 
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The Lij relate the induced charge and heat current den- 
sities to their sources such as a temperature gradient 
J9| . In the absence of interactions and inelastic scatter- 
ing processes, the Ly are expressed as jl4|,[15|[lG| 



(_l)<+i 



A(E) [E - fi(T)} l+J - 2 
df[E,n,T)- 



dE 



dE 
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for i,j — 1,2, where /x is the chemical potential of the 
system, f(E,fi,T) is the Fermi distribution function, 
and A(E) describes the system dependent features. In 
the Anderson model, one sets A(E) to be equal to the 
critical behavior of a(E) oc |1 - EjE c \ v @). Note, how- 
ever, that this behavior near the MIT does not contain 
a T dependence. Hence, the T dependence of the L^ in 
Eq. (El) is merely due to the broadening of / and the 
T dependence of \x. The latter stems from the structure 
of the density of states, variations in which yield only 
negligible changes in L^ |l3f. Thus, in order to model a 
correct T dependence of the L^ and, consequently, the 
transport properties, we need to reconsider what A(E) 
should be. 

A suitable cr(E) may be obtained from appropriate 
experimental data. The recent measurements of a by 
Waffenschmidt et al. pffl in Si:P at the MIT under uniax- 
ial stress show that a(t, T) obeys scaling with v = 1±0.1 
and z = 2.94 ± 0.3. Obtaining (i) z w d = 3 in good 
agreement with scaling arguments |4],|5| and (ii) v which 
also agrees reasonably well with the numerical results 
for noninteracting systems |IJ,[l9],gO| makes the experi- 
mental data in Ref. |L7| an excellent empirical model for 
A(E). In those experiments, t in a(t, T) is given in terms 
of the uniaxial stress and its critical value near the MIT. 
Here, we derive a functional form of <j(E) by construct- 
ing a polynomial fit or a spline curve of the experimental 
data and setting t = |1 — E/E c \. When using this data 
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Fig. 1 Scaling of thermoelectric transport properties where 
t = |1 — Ef/E c \. The different symbols denote the relative 
positions of various values of the Fermi energy Ef with re- 
spect to the mobility edge E c . 

as input for Eq. (|3j) the difference in the order of magni- 
tude in S as compared to experiments is removed fllcfl . 
Thus, following the approach of Ref . |L3| , we now study 
whether the T dependence of <r(E,T) can be scaled as 
in Eq. (Q) and whether S, K and Lq also obey scaling. 

In Fig. Ill we show that a, S, K and Lq data for dif- 
ferent t and T parameters collapse onto scaling curves 
when plotted as a function of |1 - E F /E C \/T - 31 . For 
each figure, we clearly obtain two branches, one for the 
metallic regime and another for the insulating regime. 
As depicted in Fig. |]a, a satisfies Eq. (|). With v = 1 in 
accordance with the experiment in Ref. [17]], \jvz = 0.31 
gives z = 3.2. This is in good agreement with the predic- 
tion z = 3 for 3D noninteracting systems. Furthermore, 
it indicates that the Harris criterion pi[ , vz > 1, is sat- 
isfied which in turn implies a sharp MIT. 

The prefactor in K in Fig. [j]c verifies that a and K/T 
behave similarly as the MIT is approached. This con- 
firms earlier results in the 3D Anderson model from var- 
ious methods [p[ |22[| . Meanwhile, the prefactors in both 
S and Lq in Figs. Lib and Up have not been observed in 
the respective scaling curves in the 3D Anderson model 
0. The prefactors serve as corrections to the scaling of 



S and Lq when an appropriate cr(E, T) is used as input 
in Eq. §. 

As shown in Fig. [yd in accordance also with the re- 
sults in Ref. |l3| , Lq — > 7r 2 /3, as the MIT is approached 
from the metallic or the insulating regime. This is the 
expected value in the Sommerfeld free electron theory. 
It is different from the result for the unmodified 3D An- 
derson model B for which the magnitude of Lq depends 

In conclusion, we find that by modifying a(E,T) in 
our calculations, the thermoelectric transport properties 
near the MIT obey scaling. 
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